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Abstract: 

In order to understand the initial partonic state in proton-nucleus and electron-nucleus 
collisions, we investigate the total, inelastic, and (quasi-)elastic cross sections in pA and 
7 *A collisions, as these observables are insensitive to possible collective effects in the final 
state interactions. We used as a tool the DIPSY dipole model, which is based on BFKL 
dynamics including non-leading effects, saturation, and colour interference, which we have 
extended to describe collisions of protons and virtual photons with nuclei. We present re¬ 
sults for collisions with O, Cu, and Pb nuclei, and reproduce preliminary data on the pPb 
inelastic cross section at LHC by CMS and LHCb. The large NN cross section results in pA 
scattering that scales approximately with the area. The results are compared with conven¬ 
tional Glauber model calculations, and we note that the more subtle dynamical effects are 
more easily studied in the ratios between the total, inelastic and (quasi-)elastic cross sec¬ 
tions. The smaller photon interaction makes the 7 *A collisions more closely proportional to 
A, and we see here that future electron-ion colliders would be valuable complements to the 
pA collisions in studies of dynamical effects from correlations, coherence and fluctuations 
in the initial state in high energy nuclear collisions. 
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1. Introduction 

An important question in the understanding of the strong force is the behaviour of a hot 
and dense plasma, and properties of the QCD phase diagram, with a possible critical 
point. We note that it is here very difficult to get definitive answers from ab initio lattice 
simulations [1,2]. Results from RHIC, summarised in the “White Papers” by the four 
collaborations [3-6], have been interpreted as indicating the dynamics of a nearly perfect 
fluid. These results were subsequently extended to higher energies and investigated in a 
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broader kinematic range with improved precision by the LHC experiments that have heavy 
ion programmes. These experiments include ALICE, ATLAS and CMS, and their status 
during 2014 were summarised for example by references [7-9]. The pPb program at the 
LHC turned out to be very surprising, culminating in observation of signals of collectivity 
not only in PbPb, but also in pPb reactions. Recent results from the Beam Energy Scan 
program at RHIC [10] are adding important elements to the picture in the energy range 
where many observables change together. These observations have also been interpreted 
as a signal for a critical point in the theory of strong interactions, QCD, in the vicinity of 
{TE,^iB,E) = (165,95) MeV [11]. 

R has often been suggested that studies of pA or eA scattering can play an important 
role in studies of the transition between the “simpler” pp or ep collisions and the more 
complicated AA collisions. A recent review [12] on dAu and pPb collisions at RHIC and 
LHC, respectively, indicates that many of the signatures of nearly perfect fluid dynamics 
considered for AuAu or PbPb collisions appear also in dA and pA collisions at RHIC and 
LHC energies, thus suggesting that similar perfectly flowing quark matter may be created 
also in these collisions, although its volume and its lifetime must be smaller than in the 
corresponding heavy ion collisions. The similarity between pPb and PbPb collisions is also 
emphasised in the CMS overview in ref. [9]. 

We also note that even high energy pp collisions exhibit many features interpreted as 
signals for collective behaviour, such as systematics of single-particle spectra and Bose- 
Einstein correlation radii [13], the increased production of strange particles and baryons 
(not least of strange baryons) [14-16], and signatures expected from hydrodynamic flow 
like the p/tt ratio [17] and a near side ’’ridge” [18]. 

Model predictions for final states in AA or pA collisions depend strongly on the ini¬ 
tial conditions obtained from parton-parton sub-collisions. In many analyses the initial 
conditions are estimated using the Glauber model [19,20]. Based on multiple diffractive 
sub-collisions, it automatically satisfies unitarity. It was early pointed out by Gribov [21], 
that diffractive excitation of the intermediate nucleons gives a significant contribution (on 
the 10-15% level) to Glauber’s original model. Within the Good-Walker formalism [22] 
diffractive excitation can be interpreted as a result of fluctuations in the nucleon’s partonic 
substructure [23]. It was also pointed out by Bialas, Bleszyhski and Czyz [24] that fluc¬ 
tuations and correlations in the positions of the individual nucleons within a nucleus has 
an important effect. Such fluctuations also lead to so called quasi-elastic collisions, where 
the nucleus is diffractively excited. These effects of fluctuations are most easily taken into 
account in a Monte Carlo (MC) simulation, where it is also easy to account for a realistic 
nucleon distribution within the nucleus, including correlations from a hard core in the NN 
interaction (see e.g. refs. [25-27]). 

However, based on independent nucleon-nucleon collisions the Glauber model does not 
naturally include possible interference effects between partons in different nucleons, and 
it is also difficult to include uncorrelated fluctuations between the nucleons in a target 
nucleus. We also note that although Monte Carlos including diffractive excitation are 
available, see e.g. ref. [28] and further references therein, lacking good experimental data 
for inelastic diffraction, these corrections are frequently neglected, and the calculations and 
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implementations of the Glauber model are typically based on either the total NN cross 
section or the inelastic cross section. 

The high parton densities at high energy imply, however, that saturation and interfer¬ 
ence effects ought to be important. To draw firm conclusions about a collective behaviour 
in the final state interaction, it is therefore crucial to be able to separate features of the 
initial state from the dynamics of final state interactions. The aim of this paper is to inves¬ 
tigate the initial state of a colliding nucleus in greater detail, and to isolate it from effects 
of final state interactions. To this end we focus on total, elastic, and quasi-elastic cross 
sections in p^d and eA collisions, which are insensitive to effects of final state interactions. 

Going beyond the standard Monte Garlo implementations of the Glauber model, we 
are in particular interested in the following features: 

• Effects of the partonic substructure in individual nucleons within a nucleus. 

• Interference between partons located in different nucleons, but having the same 
colour. 

• Diffractive excitation including high mass diffraction and quasi-elastic scattering. 

Saturation and correlation effects are most easily treated in transverse coordinate 
space. For our analysis we will use the so-called DIPSY model [29-32], which is based 
on BFKL evolution but includes essential non-leading-log effects, conhnement, and satu¬ 
ration within the evolution. DIPSY was originally developed for 7 *p and pp collisions, but 
can be directly extended to model collisions with nuclei [33]. Besides making it easier to 
include correlations, fluctuations, and diffraction, the formulation in transverse coordinate 
space also makes the extension to model realistic nuclei fairly straight-forward. However, 
as the gluons from different nucleons can interfere, the individual nucleons are not quite 
independent, and we therefore expect non-trivial effects of shadowing. 

The outline of the paper is as follows. First, we will in section ^ briefly go through the 
theoretical foundation of the dipole formalism and its implementation in the DIPSY Monte 
Garlo. Section ^ deals with the modifications introduced in DIPSY to simulate collisions 
including heavy ions, and the way we distribute the nucleons inside a nucleus. In section ^ 
we first demonstrate that DIPSY reproduces total pp and preliminary pH cross sections 
from RHIG to LHG energies (.ys(jvAr) = 200 GeV to 8 TeV), before we describe our results 
and predictions for the total, inelastic, and (quasi-)elastic pH and 7 *H cross sections. In 
section ^ we discuss the interpretation of the results, and end the paper with a conclusion 
and outlook in section 

2. The Lund Dipole Cascade Model DIPSY 

Most event generator models of hadronic interactions at high energies are based on multi- 
parton interactions, some examples are PYTHIA [34], HERWIG [35], DTUJET [36] and 
EPOS-LHC [37]. The strongly rising parton distributions for small x, consistent with 
BFKL evolution (albeit with important non-leading effects), make the parton-parton cross 
sections very large. At higher energies, non-linear effects therefore become increasingly 
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important, and these effects are further enhanced in collisions with nuclei. The non-linear 
effects depend on the extension of the cascades in transverse space (see the pioneering 
work by Gribov, Levin, and Ryskin (GLR) [38]). While models for hadronic collisions 
generally are formulated in momentum space, analyses of saturation and non-linear effects 
are therefore more easily treated in transverse coordinate space and rapidity. 

Formalisms for interactions with high density targets often assume a homogenous tar¬ 
get with no boundary effects (which facilitates the calculation of analytic results), and 
many applications are based on the Color Glass Condensate (CGC) [39] or the BK equa¬ 
tion [40,41]. Here also effects of correlations and fluctuations are frequently neglected. 

The DIPSY model was first developed for hadronic collisions and DIS, but could also 
be directly applied to collisions with nuclei [42]. Implemented in the form of a Monte 
Carlo event generator, it can easily take into account effects of correlations, fluctuations, 
and finite nuclear geometry. When modelling pA or AB nuclear reactions with DIPSY it is 
assumed that these interactions are dominated by absorption into inelastic channels. Thus 
the imaginary part of the 5-matrix is neglected, and the total, elastic, and quasi-elastic 
cross sections are directly obtained via the optical theorem. The assumed reality of the 

5- matrix also implies that it is possible to take the Fourier transform of the amplitude in 

6- space. This method was applied to obtain the differential elastic cross section da/dt of 
pp collisions from DIPSY in ref. [42]. 


2.1 Dipole evolution in transverse coordinate space 

The DIPSY model is based on Mueller’s dipole cascade model [43-45], which is a formulation 
of leading-log BFKL evolution [46,47] in transverse coordinate space. Mueller’s model relies 
on the fact that initial-state radiation from a colour charge (in a quark or a gluon) within 
a hadron is screened at large transverse distances by an accompanying anti-charge, and 
that gluon emissions therefore can be described in terms of colour-dipole radiation. Thus 
the partonic state is described in terms of dipoles in impact-parameter space, evolved in 
rapidity when a dipole is split into two dipoles by gluon emission. The screening implies 
a suppression of large dipoles in transverse coordinate space, which is equivalent to the 
suppression of small /cj_ in the conventional BFKL evolution in momentum space. 

For a dipole with charges at the transverse points x and y, the probability per unit 
rapidity (Y) to emit a gluon at z is given by 


d'Pg ^ « ,2 jx-y)^ 

dY 2 tt {x — zY{z — yY' 


with a 


Ncas 

TT 


( 2 . 1 ) 


The emission produces two new dipoles, (a;,z) and {z,y), which can split independently 
by further gluon emissions. Repeated emissions form a cascade, with dipoles connected in 
a chain. The chain can be thought of as dipoles connected by gluons, or gluons connected 
by dipoles [48]. When two cascades collide, a dipole (£Ci,r/^) in a right-moving cascade can 
interact with a left-moving dipole (* 2 , 2 / 2 ); with probability (in Born approximation) 


P12 


^ r, f - y2f{yi - X2)^ 
4 |_ v(®i-®2)^(yi- 2 / 2 )^ 



( 2 . 2 ) 
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We note that in this leading log approximation, the probability for dipole emission in 
eq. (2.1) diverges for small dipole sizes [x — z) or (z — y), but as the interaction probability 
in eq. (|2.2|) then goes to zero, the total cross section is finite. (This is naturally related to 
colour transparency.) Although the final result is finite, the divergence causes numerical 
problems in a leading log Monte Carlo simulation [49]. This problem is removed by the 
non-leading effects included in DIPSY, which suppress very small dipoles. 

Via the optical theorem the expression in eq. (^) also equals twice the elastic dipole- 
dipole scattering amplitude. (As mentioned above we assume that the interaction is driven 
by the inelastic absorption.) Denoting the right-moving dipoles i and the left-moving j, the 
total elastic scattering amplitude is in the Born approximation given hy F = 

Multiple scatterings are taken into account in the eikonal approximation via the unitarized 
elastic amplitude (defined by the relation S = 1 — T, which makes T real) 


r = 1 - exp (-F) = 1 - exp (- ^ Pij/‘^) 


(2.3) 


To get the full elastic amplitude we also have take the average over all possible cascades. 
In impact parameter space, the differential total and elastic cross sections are then, via the 
optical theorem, given by 


da tot/d^h = 2{T), 

daei/dH = (T)2. (2.4) 

The cross section for diffractive excitation can be obtained in the Good-Walker [22] 
formalism, which is described in appendix ^ (Diffractive excitation is usually described in 
either the Good-Walker or the triple-reggeon formalism, but as demonstrated in ref. [50], for 
high mass diffraction these are just different sides of the same BFKL evolution dynamics.) 
Diffractive excitation is here determined by the fluctuations in the scattering process. A 
proton is a linear combination of all possible parton cascades, and as discussed in refs. [31, 
51], we assume that these cascades represent the diffractive eigenstates, which thus can 
come on shell in the interaction. The cross section for diffractive excitation is thus described 
by the fluctuations among the possible parton cascades in the projectile and the target. 

As described in the Appendix, the total diffractive cross section, including elastic 
scattering, is obtained by summing over all possible diffractive states, and given by (T^). 
This is the sum of elastic scattering, single diffraction of the projectile or the target, and 
double diffraction. Single diffraction of the projectile with an elastic target is similarly 
obtained by summing over all projectile states but keeping the target intact, and then 
subtract the elastic scattering (which is still given by (T)^). The result is 

dasDjd^b={{T)l)p-{T)lt. (2.5) 

The averages are here taken over all possible cascades for the projectile (P) and the target 
(r), for a fixed impact parameter b. (Note again that the amplitude T, given by eq. (^), 
is a real quantity.) 

We note that summing over projectile cascades gives production of all possible cascades 
limited by the rapidity determined by the Lorentz frame used in the analysis. This gives 
also the possibility to calculate the differential single diffractive cross section dasn / dM‘j^. 
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Double diffraction is similarly given by total diffraction minus single diffraction and 
minus elastic scattering, which gives the result 

daT)T)/d^b = {T^)p,t - {{TfT)p - {{Tfph + {T)%t- (2-6) 

2.2 Beyond leading log 

In a series of papers [29-32] a generalisation of Mueller’s model, implemented in the Monte 
Carlo event generator DIPSY, has been described in detail. Here we will only discuss the 
main points. The basic idea behind the model is to include important non-leading effects 
in the BFKL evolution, saturation effects within the evolution, and confinement. 

The full next-to-leading logarithmic corrections to BFKL have been calculated and 
have been found to be very large [52,53]. A physical interpretation of these corrections 
has been presented by Salam [54], and a dominant part is related to energy-momentum 
conservation. In the DIPSY model this is achieved by equating the emission of a gluon 
at small transverse distances with high transverse momenta for the emitted and recoiling 
gluons. Thus the conservation of energy and momentum implies a dynamic cutoff for 
very small dipoles with correspondingly high transverse momenta. This constraint has 
also important computational advantages, removing the divergence in the dipole splitting 
probability in eq. (H- Other important non-leading effects are the running coupling, 
as{p\), and the “consistency constraint” or “the energy scale terms”, which implies that 
the emissions are ordered in both the positive and negative light-cone components [55-57]. 


2.3 Nonlinear effects 


Besides these perturbative corrections, confinement effects are included via a small gluon 
mass, and non-linear saturation effects through the so-called swing mechanism, which is a 
colour suppressed correction to the cascade evolution. Mueller’s dipole evolution is derived 
in the large limit, where each colour charge is uniquely connected to an anti-charge 
within a dipole. Saturation effects are here included as a result of multiple dipole inter¬ 
actions, in the frame used for the analysis. Such multiple interactions give dipole chains 
forming loops (see ref. [30]), and are related to multiple pomeron exchange. When multi¬ 
ple interactions can be neglected, Mueller’s model gives a result, which is independent of 
the Lorentz frame used for the calculation [58]. Loops formed within the evolution are, 
however, not included. This means that when multiple interactions cannot be neglected, 
the result is no longer frame independent. (Loops within the evolution are also neglected 
e.g. in the non-linear BK equation [40,41].) 

The dipole interaction in eq. (^) is proportional to as, and thus colour suppressed 


compared to the dipole splitting in eq. (2.1). Loop formation is therefore related to the 
possibility that two dipoles have identical colours. For two dipoles with the same colour, 
we have actually a colour quadrupole, where a colour charge is effectively screened by the 
nearest anti-colour charge. Thus approximating the field by a sum of two dipoles, one 
should preferentially combine a colour charge with a nearby anti-charge. This interference 
effect is taken into account in DIPSY in an approximate way, by allowing two dipoles with 
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Figure 1: Two dipoles with the same colour, forming a colour quadrupole, may be better approx¬ 
imated by two dipoles formed when the colour charges are recoupled, as illustrated here. 


the same colour to recouple forming new dipoles, in a way favouring small dipoles. This 
mechanism, called a swing, is illustrated in figure ||; see further ref. [30]. 

Although this recipe is an approximation, the dependence on the Lorentz frame used 
for calculating the cross section is significantly reduced, as demonstrated and detailed in 
section ^.21 . 

In the simulation the swing is handled by assigning all dipoles a colour index running 
from 1 to not allowing two dipoles connected to the same gluon to have the same 
index. A pair of two dipoles, {xi,y]) and ( 032 , 1 / 2 )) having the same colour, may be better 
approximated by the combination ( 031 ,^ 2 ) s-nd ( 032 ,^ 1 ), if these dipoles are smaller. In 
the evolution the pair is allowed to “swing” back and forth between the two possible 
configurations shown in figure ||. The swing mechanism is adjusted to give the relative 
probabilities (a3i — y 2 f‘{x 2 — yi)^ : (a3i — yif‘{x 2 — 2 / 2 )^, thus favouring the configuration 
with smallest dipoles. In the evolution in rapidity, the swing is competing with the gluon 
emission in eq. (2.1), where a Sudakov-veto algorithm [59] can be used to choose which of 
the two happens hrst. 

The non-linear GLR [38] and BK [40, 41] equations are formulated in a way where 
the number of dipoles can be reduced, via a (2 —)• 1) vertex. What is calculated in these 
formalisms is actually the interaction probability. With the swing mechanism the number 
of dipoles is unchanged. However, as the dipole swing leads to smaller average dipole sizes, 
the interaction probability in eq. ( |2.2| ) is also reduced in accordance with colour trans¬ 
parency. Thus the swing represents a similar kind of physics as these evolution equations. 
An advantage with the DIPSY formalism is that it is possible to also take into account 
correlations [60], fluctuations, and effects of the nuclear size and shape. 


2.4 Final states 

The BFKL equation and Mueller’s model describe the inclusive cross section, but as seen 
from the CCFM formalism [55,61] it is possible to include softer gluon emissions to form 
exclusive final states. As described in ref. [56] in the Linked Dipole Chain model (LDC), 
these softer gluons can be added as final state radiation from an initial BFKL ladder. When 
the two evolved systems collide, some of the dipoles in the right-moving system interact 
with some in the left-moving one, and the cross sections are calculated via eqs. (H - iM)- 
The interaction enables the gluons in the interacting dipoles to come on-shell, together with 
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all parent dipoles, while non-interacting dipoles must be regarded as virtual and thus be 
reabsorbed. In a situation, where saturation is important, it is also necessary to include 
possible swings in the final state [30, 62]. In the DIPSY MC the gluons are traced in 
both momentum and coordinate space, and therefore the model also allows us to generate 
the final-state momentum distribution of gluons and, after hadronization, final hadronic 
states [32,62]. In the present paper we restrict our study to inclusive total and inelastic 
cross sections, while final states will be studied in future work on p^ and eA collisions. 

The stochastic nature of BFKL evolution also implies, that the exclusive final states 
obtained in an event generator like DIPSY, includes fluctuations and correlations which 
are essential for many observables. As an example, we note that the fluctuations in the 


scattering amplitude determines the cross section for diffractive excitation in eq. (2.4) 


2.5 The Monte Carlo event generator DIPSY 

The Lund Dipole Cascade Model has been implemented in the Monte Carlo event generator 
called DIPSY, with applications mainly to pp collisions and DIS. In a high energy collision, 
two hadrons are evolved from their respective rest frames to a Lorentz frame in which they 
collide. In its own rest system a proton is currently modelled by a simple triangle of gluons 
connected by dipoles, and the gluonic Fock state is built by successive dipole emissions 
of virtual gluons. The small-x partons are dominated by gluons, and the result for small 
X is rather insensitive to the choice of initial parton configuration, apart from its overall 
size. We also note that the proton structure functions are well reproduced by the DIPSY 
model [63]. (Quarks are important for the final state, and here valence quarks are later 
introduced by hand). 

The parameters of the MC model have been tuned to the energy dependence of the 
total and elastic pp scattering cross sections. The best values of the tune parameters that 


govern the calculations will be shown and discussed in section 4.1. Since the system at 
very large energies is dominated by gluons, the model may use the same parameters later 
on for all kinds of reactions and observables. 

Here, we show the list of the parameters whose values substantially influence the final 
results. They are described in more detail for example in ref. [42]. 


• Rmax- Non-perturbative regularisation scale, this corresponds to the maximum dipole 
size in a given simulation, above which emissions and interactions are exponentially 
suppressed [31]. Its typical value is ps 3 GeV~^. 

• Rp\ The average size of the proton at rest. Its typical value is « 3 GeV“^. 

• Wp-. The width of the Gaussian fluctuations in proton size around Rp. As in ref. [42] 
this value was set to 0.1 GeV“^. 

• ^QCD- This is the scale parameter of as, the running coupling constant of QCD, and 
its default value in the present calculations was 0.23 GeV. 

• A^: This parameter controls the swing effect, its default value in the presented cal¬ 
culations was 1. Dipoles with the same colour are allowed to swing back and forth, 




which results in an equilibrium, where the smaller dipoles have a larger weight. = 1 
has previously [32] been shown to be large enough to reach the equilibrium. 


3. Treatment of nuclei in DIPSY 


The DIPSY model, initially developed for ep and pp collisions, can be directly generalised 
to simulate collisions with nuclei. We just have to generate random positions for all the 
nucleons within a target nucleus, and let them collide with a projectile proton, a virtual 
photon, or a similarly generated projectile nucleus. In the present paper we will focus on 
high energy p^ and 'y*A scattering, and of particular interest will here be effects of colour 
interference between different nucleons in a nucleus, and effects of fluctuations, both in the 
distribution of nucleons in a nucleus, and in the parton evolution inside individual nucleons. 

High energy nuclear collisions are usually analysed within the Glauber formalism [19, 
64]. Here it is assumed that the projectile nucleon(s) travel along straight lines and undergo 
multiple diffractive sub-collisions. Inspite of its pure geometric approach, it has been quite 
successful in describing many characteristics of reactions with nuclei, and has been widely 
used in experiments at RHIC and LHC, e.g. to estimate the number of binary nucleon- 
nucleon collisions and the number of participant nucleons as a function of centrality. In 
order to illustrate some features of our results, we will below also make comparisons with 
Glauber simulations. 

Concerning fluctuations we note that there are two different origins for fluctuations i 
collisions with nuclei. The first is due to fluctuations in the position of the nucleons in a 
nucleus, while the second is related to diffractive excitation of the wounded nucleons. As 


described in section 2.1, diffractive excitation of a proton is determined by fluctuations in 


the internal proton substructure, as given by the Good-Walker formalism. The two effects. 


and their relation to the Glauber formalism, will be described in sections 3.1 and 3.2. The 


Good-Walker formalism can also be directly applied to calculate quasi-elastic scattering. 


where the nucleus is breaking up. This is also discussed in section |3.2 . 

A very interesting question is also to what extent colour charges in one nucleon can 
screen charges in other nucleons. In DIPSY this interference effect is described by the swing 


mechanism, discussed in section 2.3, and in a nucleus we also allow dipoles in different 
nucleons to swing. The result is a kind of colour reconnection, and in the analysis of 
our results we will compare the outcome with and without this inter-nucleon swing. Such 
possible interference effects between the different nucleons are not taken into account in 
the Glauber approach. 


3.1 Distribution of nucleons in a nucleus 

It was early pointed out by Bialas, Bleszyhski and Czyz [24] that fluctuations and correla¬ 
tions in the positions of the individual nucleons within a nucleus has an important effect. 
This effect, which is most easily taken into account using MG simulation techniques (not 
easily available at the time) has been studied by several authors, see e.g. refs. [27,65]. The 
nucleon positions are typically generated so as to reproduce the charge density observed in 
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eA scattering, and the hard core in the NN interaction has an important effect reducing 
the fluctuations in the distributions. 

The GLISSANDO method 

We will in our simulations follow the method in the GLISSANDO MC [25,26], to 
generate the nucleon distributions. The nucleon positions are generated according to a 
Wood-Saxon distribution for the nucleons, with the following form: 


_ pq{1 + wr"^ / B?) 

1 + exp((r — R)/a)' 


(3.1) 


Here R is the nuclear radius, a is “skin width”, and po is the central density. The parameter 
w describes a possible non-constant density, but is zero in the fits to nuclei used in this 
paper. 

The nucleon centres are randomly generated in such a way that the charge distribu¬ 
tion determined in ref. [66] is recovered, when the result is convoluted with the charge 
distribution within nucleons. The nucleons are also generated with a hard-core, which thus 
introduces short range correlations among the nucleons. As shown by Rybczynski and 
Broniowski [67], the correct two-particle correlation can be obtained if the nucleons are 
generated with a minimum distance equal to 2rcore- 

We note that if the nucleons are generated within a specific volume, the resulting 
distribution will (for a finite nucleus) be confined within a smaller volume, and its centre 
will be shifted. According to ref. [26], the correct charge distribution is, for mass numbers 
A > 16, obtained using randomly generated nucleon centres described by the Wood-Saxon 
form in eq. (3.1) with the following parameters: 


= (1.1H^/3 _ 0.656^-1/3) a = 0.459 fm, w = 0, 


(3.2) 


together with a hard core with radius Vcore = 0.45 fm. 

Currently the DIPSY MC includes parametrisations for He, O, Cu, Au, and Pb^. (Thus 
we use the same spherical form also for the light nuclei He and O.) For the nuclei studied 
in this paper this corresponds to the following radii: R^q = 6.406 fm, R^'q = 6.288 fm, 
R%'^ = 4.236 fm, and R%(j = 2.511 fm. 


3.2 Diffractive excitation of wounded nucleons 


In Glauber’s initial formulation the nucleus was described by a smooth distribution deter¬ 
mining the absorption probability, but it was early pointed out by Gribov [21], that effects 
of diffractive excitation of the wounded nucleons are quite significant. Diffractive excita¬ 
tion occurs when the projectile is a linear combination of states with different absorption 
probabilities. Diffractive excitation is here determined by the fluctuations in the scattering 


amplitude, as seen in eq. (lA.lOf) . At lower energies this could be well approximated by a 
single diffractive state, but it is now well known that diffractive excitation in pp collisions 
or DIS is not limited to low masses, and has a high cross section ~ 10 — 15% of the to¬ 
tal [68-70]. An important point is here that, due to time dilation, in a pA collision the 


^ Other nuclei can easily be added by the user. 
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state of the projectile is frozen between the sub-collisions, while the target nucleons may 
all be in different states [71]. 

As mentioned in the introduction, this effect has been studied in several analyses and, 
as an example, a model by Strikman and coworkers [28,72] has been applied by the Atlas 
collaboration [73]. However, lacking good experimental data for inelastic diffraction, these 
corrections are frequently neglected, with calculations based on either the total NN cross 
section or the inelastic cross section (including or excluding diffraction). 

In the simplest form of the Glauber Model, the target in a NN collision acts as a black 
absorber. The projectile nucleon travels along a straight line and interacts inelastically if 
the transverse distance to a nucleon in the target is smaller than a distance R, with 


ttR^ 


= a 


NN 

in,ND' 


(3.3) 


Here 

disc approximation 
eqs. (|2^, pH^)): 


inelastic, non-diffractive nucleon-nucleon cross section. This black- 
implies the following cross sections for nucleon-nucleon collisions (see 


cTtot = 2 J <fb{T{h)) = 27rR‘^ 

(Tel = J d?h{T{b))‘^ = -kR? 

an = I d%mbf) - {T{b)f) = 0 
fTin,ND = j d^{l - (1 - T(6))2) = 7rR\ (3.4) 


These relations demonstrate, that in a black-disc approximation the diffractive part of the 
cross section vanishes, and that the elastic and inelastic cross sections are equal, both given 
by half the total. 

The most simple model accounting for diffractive excitation is the so called grey disc 
model. Here it is assumed that within a radius R the projectile is absorbed with probability 
a, with 0 < a < 1. The resulting pp cross sections are here 

= 2 J d^b{T{h)) = 27TR^a 

= j d%{T{b)f = TTR^a^ 

= J dH{{T{hf) - {T{b)f) = irR^il - a) 

= J d%{l - (1 - T{b)f) = nR^a. (3.5) 

The parameters R and a can now be adjusted to reproduce e.g. the total and the elastic pp 
cross sections. It is interesting to note that even in this gray disc limit, crei + c’'D = <Tin,ND = 
(Ttot/2. We also note that in the LHC energy region ~ 0.25(Ttot) which in this grey disc 
model gives (Td ~ (Tei and crin,ND ~ 0.5(Ttot- Also other forms for the pp interaction have 
been considered in the literature, e.g. a Gaussian interaction profile [64,74,75]. 


O’tot 

<Tel 

<TD 

<Tin,ND 
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In the black disc approximation the fluctuations are totally neglected. As mentioned 
above, in a pA collision the state of the projectile is frozen between the sub-collisions, while 
the target nucleons may all be in different states. Averaging over different colliding target 
nucleons therefore tends to reduce the fluctuations. Consequently, if the parameters in the 
grey disc approximation are adjusted to the total and elastic pp cross sections, the effects 
of fluctuations are instead overestimated. The averaging over different states for the target 
nucleons is properly taken into account in DIPSY, and when comparing the DIPSY results 
with Glauber calculations we therefore expect the DIPSY to lie in between the black and 
grey disc results 

As the main difference in pA collisions between the DIPSY model and Glauber-based 
models is the treatment of fluctuations inside nucleons, it is interesting to study the diffrac¬ 
tive dissociation of the nucleus. From the single diffraction cross section in eq. ( |2.5D we 
immediately get that this is given by 

dasB,A/d^b = {{T)1 )a - {T)Ia ■ (3-6) 

Experimentally this means that we must detect an elastically scattered proton in one 
direction, and a dissociated nucleus in the other. In DIPSY there is in this case a dependence 
on the frame in which the collision is studied. Typically the proton-nucleon rest frame is 
used, which means that we only consider the case where the diffractively dissociated nucleus 
does not produce particles in the proton hemisphere of the detector. 

It should be noted that even if the grey- and black-disc Glauber calculations do not 
model the single diffractive components for pp collisions, they will give diffr actively disso¬ 
ciated nuclei due to the fluctuations in the distribution of nucleons. This can be thought 
of as elastic scattering of the proton with one of the nucleons, causing the nucleus to break 
up. 

Experimentally it may be difficult to measure low-mass diffraction in the nucleus di¬ 
rection, and below we will also study what we call quasi-elastic scattering, where again the 
proton is scattered elastically, but the only other requirement is that the proton hemisphere 
of the detector is empty. This corresponds to adding the single diffractive and elastic cross 
sections, resulting in 

da,u/dH = ml)A ■ (3.7) 


4. Monte Carlo cross section results 

In this section we present various Monte Garlo simulation results for high energy proton- 
proton, proton-ion, photon-proton and photon-ion reaction cross sections performed by 
DIPSY. In order to show that these numbers are reasonable, they are compared to (prelim¬ 
inary) data, where these data are already available. Most of the (heavy) ion results (pA 
and 7 *A) are in fact predictions, given that the corresponding experimental values are not 
yet measured and many of these determinations are not even foreseen at present. 

^This overestimate of the fluctuations in the target is also present in other Glauber approaches, where 
the same diffractive eigenstate is used in each NN sub-collision. 
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As discussed in section |l|, our calculations are based on the evolution and interaction of 
virtual gluons whose number is increasing with the increase of ^/s and the size of the nuclei 
involved. As a practical consequence, the simulations for more energetic reactions and 
higher values of A are limited by the available computational power to achieve reasonable 
statistics in reasonable time. The scope of this paper, due to this practical consideration, 
covers ranges of centre-of-mass energies from the upper range of RHIC energies of y/sjsfjsf = 
200 GeV up to the LHC energies of 8 TeV, selecting a few typical points in between. The 
RHIC energy ^snn = 200 GeV for colliding nuclei provides partons with about x = 0.01, 
which is the lowest limit for an acceptable precision with the event generator. At lower 
energies the contribution of the previously neglected quark degrees of freedom may play 
quite a significant role. The upper range of the energy scale was chosen to be 8 TeV so that 
the simulations could be completed with the available computing power in the timescale of 
this project. For the particular case of the pPb reactions, we have evaluated the total cross 
section and its decomposition also for the expected, future LHC energy of y'sNN = 10 TeV. 

In the figures, as a general rule, DIPSY results are shown for eight different nucleon- 
nucleon centre-of-mass energies chosen appropriately from RHIC to LHC energies in case 
of pp and pA reactions. To help the visualisation of the results, the simulated points 
were connected by simple straight lines to guide the eye. For 7 *Asimulations, one more 
simulation point was added at y/TjTv = 100 GeV for the sake of possible eA experiments 
at RHIC where, in the current planning phase, the = 20 — 90 GeV centre of mass energy 
range is considered for eA collisions with large A, and the 30 to 145 GeV energy range is 
considered for polarised ep collisions [76,77]. 

4.1 Tuning and comparison to data 

4.1.1 pp total and elastic cross sections 

In figure DIPSY results of pp total and elastic cross sections are shown after the most 
sensitive model parameters are tuned to data. The ball park values of these parameters 
are listed in section |2.5| , and we give their actual values in this sub-section, with the help 
of these tunes. Subsequently, these parameters were fixed and utilised in the more complex 
simulations of proton-ion and photon-ion events. For comparisons they were also used to 
simulate photon-proton reactions. 

Our plot in figure follows the style of figure 1 of ref. [78], where the TOTEM col¬ 
laboration compared its recently measured total and elastic pp cross sections at 7 and 
8 TeV with earlier results, and low statistics cosmic ray experiments at very high energies, 
see refs. [78, 79] for further details. Cross sections measured at lower energies were taken 
from from the Particle Data Group [80] database. The best fits with a formula by the 
COMPETE Collaboration [81] are also compared to the data on the total and also on the 
elastic cross sections with coefficients shown for cJei in the legend of the plot. The results of 
the simulations by DIPSY (shown by full dots) indicate that DIPSY simulations follow the 
trends of the cross-section data reasonably well. 

To achieve this reasonably good level of description of the total and elastic cross sec¬ 
tions, the QCD scale parameter Aqcd in DIPSY was tuned to the value of 0.230 MeV. The 
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pp and pp cross sections 



Figure 2: Total and elastic cross sections in pp and pp collisions plotted together with tuned DIPSY 
simulation results (red filled circles). The total cross section data (open triangles) are described by 
the best fit of the COMPETE collaboration [81] and the inelastic cross sections are also fitted by 
a phenomenological formula as shown in the figure. 


effective size parameter of the proton was found to be Rp = 2.9 GeV ^ ~ 0.57 fm (which 
is different from 0.45 fm used for the hard-core radius mentioned above). 

The successful reproduction of the experimental data on the total and elastic cross 
sections of pp reactions in the energy range of y/s = 0.2 — 8 TeV energy range, where DIPSY 
has valid approximate assumptions for the parton evolution within the nucleons, opens the 
way for Monte Carlo calculations of electron-ion (eA) and proton-ion (pA) reactions. 
Heavy ion {AB) collisions can be studied in a similar manner, but the investigation of 
these reactions goes beyond the scope of this article. 


4.1.2 pPb inelastic cross sections 


Above the highest RHIC ion collision energy of y/SMN = 200 GeV, in the region where 
DIPSY can be reliably applied, there are only a few heavy ion cross sections measured. One 
of them is the pPb total inelastic cross section for which the CMS [82] and LHCb [83] 
collaborations have presented preliminary results^ as shown in figure The DIPSY results, 
shown by a solid line connecting the simulated points, agrees well with the data, as does 
the two additional theoretical calculations shown from the EPOS-LHC [85] and QGSJetll- 
04 [86] models. 

The DIPSY results are mainly driven by the underlying geometrical description of the 


nucleus from the GLISSANDO parameterisation, as discussed in more detail in section 4.3 


®Note that the LHCb value required events with at least one track in their acceptance region, while the 
CMS value has been fully corrected to the total inelastic cross section. Also ALICE has presented a related 
result for visible cross section [84] consistent with LHCb. 


- 14 - 











2.5 


pPb inelastic cross sections 



Figure 3: pPb inelastic cross sections recently measured by the CMS (black triangle) and the 
LHCb (blue square) collaborations. The solid line guides the eye for the DIPSY results starting 
from RHIC collision energies. The predictions of QGSJetII-04 [86] and EPOS-LHC [85] simulations 
are indicated by dotted and dashed lines, respectively. 

4.2 Predictions for pA and 7 *A total cross sections 

4.2.1 pA total cross sections 

In figure § total cross sections calculated by DIPSY are shown for collisions of protons 
with different types of ions. Lead, copper and oxygen were selected to represent heavy, 
intermediate, and light ions to get an overview about the main characteristics of pA cross 
sections. We note that these cross sections grow more slowly with energy than the pp cross 
section in figure p. This is more clearly seen in figure |^, where we show the ratio of the 
total pA cross section to the corresponding total pp cross section. For clarity, we here also 
normalise this ratio by the mass number, A, for the colliding ion. If the NN cross section 
would be small, this ration should be one. The deviation from one, which is larger for the 
heavy lead nucleus and for higher energy, is thus a related to the average number of 
collisions in a single pA event. 

For a completely black nucleus, the cross section is expected to scale with the area 
overlap for the two projectiles, i.e. oc {Ra + Rp)^■ Approximating Ra A^/^Rp will thus 
lead to a scaling ~ i?p(A^/^ + 1)^. Figure § also shows the pA total cross sections scaled 
by (A^/^ + 1)^. We can here see that, as expected, this area scaling works better for heavy 
nuclei and high energy, when the black limit is approached. Analogously the contribution 
from volume effects oc A are more visible for smaller nuclei and lower energy. 

4.2.2 7 *A total cross sections 

When a (virtual) photon fluctuates into a qq pair, it forms a colour dipole, which interacts 
with the nucleus in the same way as a dipole in a colliding proton. Thus it produces a 
dipole cascade, which is comparatively long-lived, and frozen during the interaction with 
the nucleus. In figures we show the total cross section for 7 *A collisions in the energy 
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pA total cross sections 


+ pPb 
+ pCu 
4 - i pO 



^[TeV] 


Figure 4: Total cross sections for collisions obtained from the DIPSY Monte Carlo model 
calculations. The results are shown as a function of the nucleon-nucleon centre-of-mass energy, 
covering a range from the upper RHIC energies to the available LHC energies. 



Figure 5: DIPSY predictions for the total pA cross sections scaled by (left)j and by (H^/^ + 1)^ 

(right). For the mass-number scaling the ratio for pO collisions closer to 1 than for pPb collisions, 
indicating, that volume effects are more important for lighter nuclei and at lower colliding energies, 
as far as the non-linear nuclear effects implemented in DIPSY are considered, while the opposite is 
found for the surface scaling. 


range between 0.1 and 7 TeV, similar to the energy range studied above for pA 

collisions. Results are presented for photon virtualities = 2.5, 5, 25 and 100 GeV^ 
colliding with O, Cu, and Au nuclei. We also present the results after scaling with the 7 *p 
cross section and the mass number A. We have used the parameters tuned to the pp cross 
sections in section O, and we have checked that we reproduce data on 7 *p cross sections 
to a reasonable degree. This means, however, that the 7 *A cross sections presented here 
should not be thought of as true precision predictions, and we will here only discuss the 
qualitative behavior of the cross sections and their ratios. 

We note that the dipole-nucleon cross section is suppressed compared to the NN 
interaction, as an effect of colour transparency. Therefore the cross sections scale more 
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y*0 total cross sections 




^[TeV] 


Figure 6: Left panel indicates DIPSY predictions for the total cross sections of 7*0 reactions as 
a function of the photon virtuality and the available energy n- The right panel compares 
these calculations to 7 *p reactions, scaled also by the mass number A. For comparison we also 
show the ratio pO/pp. 


y*Cu total cross sections 



y*Cu/y*p totai cross section ratios 



Figure 7: Same as 


figure H but for Cu nuclei. 


closely with the volume ~ A than with the area, as was the case for pA collisions. For 
the same reason it also grows with energy, approximately proportional to the 7 *p cross 
section. As seen in the figures, this scaling is most clear for larger and smaller A, while 
saturation effects become more noticeable for lower and heavier nuclei. 

The 7 *A cross section is also suppressed by the factor oem in the photon-quark 
coupling. Note that the cross sections are here measured in millibarn, instead of barn, 
which was used above for the pA cross sections. However, as the dipole state is frozen 
during the collision, this extra suppression has no effect on the scaling behaviour discussed 
here. 

Figure ^ indicates the results of a DIPSY simulation for the intermediately heavy Cu 
nuclei that can be compared in a straightforward way with the results shown on figure |^. 
From this comparison one can observe, that /Aa^o^ ■ This suggests 

that even for intermediately heavy nuclei such as Cu, it is more difficult for the dipoles to 
find a hole and to pass through these nuclei, as compared to lighter nuclei such as O. 
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y*Au total cross sections 




^[TeV] 

Figure 8: Same as figure ^ but for Au nuclei. 

Figure indicates that in 7 *Au reactions the mass number or volume scaling is even 
worse approximation than in 7 *Cu reactions. This indicates that for large nuclei and for 
large ^/s the nuclei become asymptotically black even for a small dipole coming from a 
virtual photon, 7 *. Nevertheless, when compared to calculations with similar centre-of- 
mass energies for pPb collisions, one can observe that virtual photons can pass through 
these large nuclei Au or Pb more easily as compared to the collision of protons with the 
same target. 

4.3 Semi-inclusive pA cross sections and comparisons with Glauber models 

Besides colour coherence effects, to be discussed in section the most important feature 
of the DIPSY model for nucleus collisions is its treatment of correlations, and fluctuations, 
not only between the nucleons, but also inside the nucleons and between partons in different 
nucleons. To investigate the observable consequences of this we will in this section compare 
some results for semi-inclusive cross sections in DIPSY with those of standard Glauber 
calculations, where mainly fluctuations in the positions of nucleons and short range nuclear 
correlations between nucleons are considered. 

In figure ^ we present DIPSY predictions for the inelastic, non-diffractive cross sections 
for pA collisions as a function of collision energy, comparing to a Glauber calculation using 
the black disc with the same nucleon distribution (see section P) as input. Glearly the 
results are quite similar with differences only at the percent level. However, if we look at 
total pA cross sections in figure the differences are much larger, almost up to 30% for 
small nuclei and lower energies. The differences here are mainly in the way how fluctuations 
are treated. 

It should be noted that the input radius to the Glauber calculation is here taken from 
the inelastic, non-diffractive pp cross sections, c.f. eq. (3.3), using the value obtained from 
DIPSY. However, looking at eq. (3.4), it is clear that within the black disc approximation 
we could also have taken the total or elastic pp cross sections to define the input radius. 
In fact, if we had taken = 27rR‘^ to define the input radius, the difference w.r.t. DIPSY 
would be much smaller. On the other hand the differences for the inelastic, non-diffractive 
cross sections would then be larger. 
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Figure 9: On the left panel Glauber Monte-Carlo and DIPSY calculations are compared for the 
^in^ND inelastic non-diffractive cross sections. In these calculations, the DIPSY results for crin,ND 
pp cross sections were used, as input, in the Glauber calculations. The Glauber Model MG results 
were obtained in the black-disc approach. The right panel shows the ratios of the DIPSY calculation 
to the Glauber Monte Carlo, indicating that the effects of correlations and fluctuations contribute 
to these variables on the typical 1% scale. 
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Figure 10: Left panel indicates a comparison of total cross section calculations from DIPSY 
to the corresponding Glauber Model MG results, using the black-disc approach. For Glauber, the 
same (Tm.ND PP cross sections were used, as input, that DIPSY calculated. Right panel indicates 
the ratios of total cross sections calculated by DIPSY to the Glauber Model MG in black-disc 
approach, indicating that although the inelastic cross sections are within 1% the same in the two 
approach, the predictions for the total cross section differ by a rather significant value that is about 
20 — 30% at RHIG energies and decreases to 10 — 15% at LHC energies. 


As pointed out in section |3.2| , the observables most sensitive to fluctuations are related 
to diffraction. For the black-disc Glauber calculations in eq. (|3.4D , the pp diffractive cross 
section is zero. For pA collisions, however, the Glauber calculation will result in a diffractive 
component, which is then entirely related to the fluctuations in the nucleon distribution 
within the nucleus. This can also be thought of as the elastic scattering of the proton 
projectile with one or more of the nucleons, resulting in a diffractive excitation of the 
nucleus. It is therefore interesting to compare this to DIPSY, where more fluctuations are 
taken into account. 

In table || we present predictions from DIPSY for different (semi-) inclusive pPb cross 
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Model 


DIPSY 

Black disc 

Black disc 

Black disc 

Grey disc 





(o-tot) 


.n) 

(o-in 

nd) 

(*^tot; ^el) 

yJSNN (T'eV) 

5 

10 

5 

10 

5 

10 

5 

10 

5 

10 

(^tot 

(b) 

3.54 

3.62 

3.50 

3.58 

3.88 

3.95 

3.73 

3.80 

3.69 

3.77 

f^in 

(b) 

2.04 

2.07 

1.95 

1.98 

2.14 

2.17 

2.06 

2.09 

2.07 

2.11 

0'in,ND 

(b) 

1.89 

1.92 

1.75 

1.79 

1.94 

1.98 

1.86 

1.90 

1.84 

1.89 

^el 

(b) 

1.51 

1.55 

1.55 

1.60 

1.73 

1.78 

1.66 

1.70 

1.62 

1.66 

O'SD.A 

(b) 

0.085 

0.086 

0.198 

0.192 

0.204 

0.198 

0.200 

0.195 

0.083 

0.085 

O'SD.p 

(b) 

0.023 

0.024 

- 

- 

- 

- 

- 

- 

- 

- 

C^DD 

(b) 

0.038 

0.038 

- 

- 

- 

- 

- 

- 

0.142 

0.137 

^el=k 

(b) 

1.59 

1.64 

1.75 

1.79 

1.94 

1.98 

1.86 

1.90 

1.70 

1.75 

^el*/^in 


0.78 

0.79 

0.90 

0.90 

0.91 

0.91 

0.90 

0.91 

0.82 

0.83 

0'in,ND/<7tot 

0.53 

0.53 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 


Table 1: DIPSY predictions for different (semi-) inclusive pPb cross sections (given in units of barns) 
at LHC energies = 5 and 10 TeV) compared to different Glauber calculations. The input 

radius to the black disc calculations were given by ttR"^ = crP^No- and crPP. For the Grey disc 

both and were used (see eq. (3.5)). In all cases using these pp cross sections obtained from 
DIPSY. The cross sections shown are the total (utot), the inelastic including diffractive excitation 
(cTin), the inelastic non-diffractive ((Tin,ND)i the elastic (crei)i the single diffractive excitation of the 
nucleus (ctsd.a) and proton (ctsd.p)) the double diffractive excitation (udd)) and the quasi-elastic 
scattering (aei* in eq. &)■ Also shown are the cross section ratios, aei*/(Jin and crin,ND/CTtot- The 
statistical uncertainty in the calculations are below a percent, except for the diffractive excitation 
cross sections, where the uncertainty is around 3%. 


sections at LHC energies (V« 5 and 10 TeV) compared to different Glauber calcula¬ 
tions. The input for the black disc radii were taken from the pp cross sections calculated 
by DIPSY, and for all cross sections we find a trivial increase as the radius becomes larger 

{aZ/2 < C; 


in,ND 


< As noted before the inelastic, non-diffractive cross section is 

close to DIPSY, when using the pp inelastic, non-diffractive cross section for the black-disc 
radius. We note, however, that the predicted cross sections for the diffractive dissociation 
of the nucleus are quite insensitive to the input radius, and are for a black disc in all three 
cases much larger than the predictions from DIPSY. This effect also shows up in the quasi¬ 
elastic cross sections (c./. section 3.2). The reason for this can be found in the fact that all 
black-disc calculations over-estimate the elastic nucleon-nucleon cross section, resulting in 
a larger probability to quasi-elastically excite the nucleus. 

It is therefore interesting to also compare with a grey-disc Glauber calculation 
{c.f. eq. (3.5)), where both the elastic and the total NN cross sections are set to the 
same as in DIPSY. In table ^ we see that the nucleus dissociation cross section in this case 
is very close to DIPSY. However, since the elastic pPb cross section comes out slightly larger 
than in DIPSY, also the quasi-elastic cross section is larger. 

The inelastic pPb cross section has already been measured at LHC [82,83], and there 
is also a possibility to measure the quasi-elastic cross section using e.g. the TOTEM ex¬ 
periment. In table |l] we also present predictions for the ratio between these cross sections. 
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as one may expect that some systematic uncertainties in the experimental measurements 
{e.g. uncertainties in luminosity) may cancel in such a ratio. 

For completeness we also show in table [I| the DIPSY results for single diffractive excita¬ 
tion of the proton, which is absent in all the Glauber calculations, and the double diffraction 
cross section, which is zero in the black disc approximations. Finally, we note that the ratio 
o'in.ND/o'tot by construction is one half in the Glauber calculations (see eqs. ( p.4D and 
while in DIPSY, the ratio is slightly higher. 

5. Discussion 

5.1 Saturation effects 

Colour reconnection (swing) is one of the nuclear effects introduced in the Lund Dipole 
Cascade Model DIPSY. Such a swing may take place if dipoles of the same colour are close 
enough in a longer chain of dipoles, or, if the colour conhguration of two different chains 
allows for a colour reconnection. This reconnection may happen with certain probability 
discussed in section Q in such a way that preferably smaller dipole chains are formed, possi¬ 
bly giving rise to chains forming small closed loops. Smaller dipoles reduce the interaction 
probability and hence result in smaller cross sections. This feature effectively corresponds 
to a gluon saturation effect, as mentioned already at the description of the DIPSY model 
in section |2|. 

In this section, we quantify the magnitude of the swing effect for electron-nucleus 
(more precisely virtual photon-nucleus) and proton-nucleus collisions. In particular, we 
switch on and off the inter-nucleon swings. However, it is important to note that we keep 
the intro-nucleon swings turned on so that the values of the total and elastic pp cross 
sections remain the same as indicated in figure 

It is well known that Gribov corrections to multiple diffractive scattering [21] are 
typically on the 10 — 15% scale and in principle they take into account the fluctuating size 
of the projectile when evaluating, for example, the total cross section of proton-nucleus or 
deuteron-nucleus collisions. In this picture, it is important to consider that the projectile 
itself has a fluctuating size, and when it is in a small size conhguration, it can pass through 
the nucleus easier. Hence the total pH cross section is typically reduced when the Gribov 
correction is taken into account. Intra-nucleon swings may increase the fluctuations in the 
size of the nucleon, and are hence related to these Gribov corrections. The inter-nucleon 
swing mechanism is an additional, new effect, that may decrease the probability of the 
interaction further, reducing the impact of large colour dipole conhgurations even further. 

In hgure |Tl| the total cross sections with and without inter-nucleon swing are compared 
for 7*0, 7 *Au, pO and pPb reactions. As the number of possible swings will increase with 
the number of dipoles, and hence with the collision energy, one might assume that the 
effects of the inter-nucleon swing would also increase with energy, however from hgure 
it is clear that the energy dependence is practically hat. From this we can conclude that 
the main effect comes from reconnections related to the initial dipoles of the nucleons. For 
a virtual photon probe we hnd the largest effect for a heavy nucleus (over 10% for 7 *Au), 
which is natural, since here we have more nucleons overlapping in impact-parameter space. 


(3.5)), 
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Total cross section ratios 



Figure 11: Total cross section ratios for pPb and pO reactions, as well as for 7*0 and 7 *An, with 
inter-nucleon of swing mechanism switched off as compared to on. 


For a proton projectile, however, the effects overall as well as the differences between nuclei 
are smaller, since for central collisions (where the overlap is largest) at these energies, the 
nucleus is already effectively “black”, while the same nucleus is much more transparent for 
a virtual photon probe. 


5.2 Frame (in)dependence as a function of 

In this sub-section we demonstrate how can one minimise the dependence of the results of 
the DIPSY model on the frame chosen as reference frame for the simulation. The reason 
for this discussion is that the Lund Cascade Model DIPSY is not formulated in an entirely 
boost-invariant manner, yet, so the simulation results depend slightly on the Lorentz-frame 
chosen as a reference frame. 

The results presented in section ^ are obtained using the NN (or 'y*N) centre-of-mass 
frame, with parameters tuned to pp total and elastic cross sections in this frame. The 
screening effects in interactions with nuclei are most sensitive to the cross section for the 
individual NN collisions. Therefore, although the result depends on the frame used when 
the parameter values are fixed, as is shown in the left panel of figure 1^, the result is 
insensitive to the frame if the parameters are tuned, in the same frame, to pp cross section 
data. This is illustrated on the right panel of figure |l^, which shows as function 

of the pp total cross section. 

This way the error due to the frame dependence is reduced to a level of about 1% 
or less, as illustrated in the right panel of figure IJ, which negligible as compared to the 
10 — 15% size of the dynamical swing effect that we explored in the previous sub-section. 


6. Conclusions 

Analyses of collective behaviour from final state properties in AA collisions depend crucially 
on the assumptions about the initial state, obtained from parton-parton interactions. This 


- 22 - 












0.5 


pPb/pp total cross section ratios 


pPb/pp total cross section ratios 



\i5;;;;[Tev] 


ofol [mb] 


Figure 12: Left panel indicates the frame dependence of the normalised total cross 

section ratios calculated in three different frames,given by yframe = — Y) ■ Vp + Y ■ yA, with 

Y = 0.5, 0.7, and 0.9, as a function of the centre-of-mass energy. We see that even these normalised 
ratios have a frame dependence of about 25%. The right panel, however, illustrates that most of 
the frame dependence enters through the frame dependence of the total pp cross section: if these 
values are correctly tuned to the measurements in the frame of the calculation, the results (with 
less than 1% error) become approximately frame independent. 


initial state is sensitive to all possible correlations, fluctuations, and colour interference 
effects between the partons within a nucleon as well as in different nucleons. To draw 
stable conclusions, it is therefore important to have a good understanding of this initial 
state. It has often been suggested that studies of p^ collisions, as intermediate between 
AA and pp collisions, will be of great help here. 

In this paper we will investigate the properties of a high energy nucleus via studies of 
total, elastic, and quasi-elastic cross sections in p^ and 7 *^ collisions, as these observables 
depend only upon the initial states, but are insensitive to any collective effects in the final 
state interactions. The Lund Dipole Cascade model DIPSY offers special possibilities to 
study the evolution of gluons inside hadrons and nuclei, including effects of unitarization 
at small x, interference, correlations, and fluctuations. The model is based on BFKL 
evolution, including essential non-leading effects, saturation, and colour interference effects. 
It reproduces successfully total, elastic, and diffractive cross sections in pp collisions and 
DIS, within the energy range from RHIC to LHC for pp collisions and the HERA range 
for DIS, and it can be directly generalised to study collisions with nuclei. 

The general features of the results are, as expected, that in pA interactions the centre is 
rather black, which means that the cross section approximately scales with the area {Ra + 
Rp)^ ~ -b 1)^. In 7 *A collisions the nucleus is much more transparent, in particular 

for high and lighter nuclei. The cross sections therefore here scale approximately with 
A, but behaves somewhat more similar to pA for smaller and/or larger A and higher 
energy. As a particularly interesting effect we note that colour coherence effects between 
gluons in different nucleons in the nucleus, can reduce the 7 *Au cross section by about 
10%. These results underline the importance of future electron-ion collider experiments to 
study the initial partonic state of highly energetic nuclei. 

To study the more subtle effects related to colour coherence and fluctuations, we make 
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comparisons with Glauber calculations. We here see that the interesting effects show up in 
the relative size of the total, elastic, diffractive, and quasi-elastic cross sections. We point 
out that both the experimental and the theoretical accuracy will be higher for these ratios 
than for the individual cross sections. 

Finally we note that the DIPSY program is also able to produce exclusive hadronic 
final states in pp collisions [32]. The extension of this to also model final states of collisions 
involving heavy ions is certainly feasible, but the effects of the extremely crowded partonic 
environment in pA and AA are not trivial and must be modelled with great care. Progress 
towards understanding string fragmentation in a dense environment has already been made 
in [87], but the full picture of the final states in the DIPSY dipole model is yet to come. 

A. Appendix. The Good—Walker formalism for diffractive excitation. 

Diffraction is normally interpreted as elastic scattering driven by absorption. If the absorp¬ 
tion probability in Born approximation is given by 2F, then rescattering exponentiates in 
impact parameter space, giving the unitarized absorption probability 

Pabs = dain/cfb = I - e~‘^^. (A.l) 

We define the elastic amplitude, T, via the relation S = 1 — T, which makes T real. The 
optical theorem, T = (l/2)(r^ -|- Pabs)^ then gives the following result for the amplitude 
and the elastic and total cross sections: 


T = l-e~^ (A.2) 

daei/(fb = (1 — (A.3) 

datotld'^b = 2(1 - e~^). (A.4) 

For a projectile with a substructure, the mass eigenstates can differ from the eigenstates 
of diffraction, i.e. states with a specific absorption probability [22]. Call the diffractive 
eigenstates 4>„, with elastic scattering amplitudes T^. The mass eigenstates are linear 
combinations of the states <I>„; 

'^k = '^Ckn^n (with = 4'i). (A.5) 

n 

The elastic scattering amplitude is given by 

(4/i|r|4/i) = ^cLTn = (r), (A.6) 

and the elastic cross section 

da,,/dH={Y,clj'n? = {T)\ (A.7) 

The amplitude for diffractive transition to the mass eigenstate is given by 

(4'fc|r|4'i) = ^ CknTnCln, (A.8) 
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which gives a total diffractive cross section (including elastic scattering) 

dcTdis/dH = = (r2). (A.9) 

k 

Consequently the cross section for diffractive excitation is given by the fluctuations: 

daji/dH = daAm/d% - da^/dH = (T^) - {Tf. (A.IO) 

If also the target has a substructure, it is possible to have either single excitation of 
the projectile, of the target, or double diffractive excitation, with cross sections given in 
section p.l| . 
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